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Abstract. 

Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations shorter and longer 
than about 2 seconds. There are, however, some indications for the existence of a third type of burst. Swift satellite detectors 
have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore it is worth to reanalyze the durations 
and their distribution and also the classification of GRBs. Using The First BAT Catalog the maximum likelihood estimation 
was used to analyzed the duration distribution of GRBs. The three log-normal fit is significantly (99.54% probability) better 
than the two for the duration distribution. Monte-Carlo simulations also confirm this probability (99.2%). 
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INTRODUCTION 

It has been a great challenge to classify gamma-ray bursts (GRBs). Using The First BATSE Catalog, [1] found a 
bimodality in the distribution of the logarithms of the durations. In that paper they used the parameter Tgo (the time in 
which 90% of the fluence is accumulated [1]) to characterize the duration of GRBs. Today it is widely accepted that 
the physics of these two groups are different, and these two kinds of GRBs are different phenomena [2, 3]. In the Swift 
database the measured redshift distribution for the two groups are also different, for short burst the median is 0.4 [4] 
and for the long ones it is 2.4 [5]. 

In a previous paper using the Third BATSE Catalog [6] have shown that the duration (Tgo) distribution of GRBs 
observed by BATSE could be well fitted by a sum of three log-normal distributions. We find it statistically unlikely 
(with a probability ^ 10^^) that there are only two groups. Simultaneously, [7] report the finding (in a multidimen- 
sional parameter space) of a very similar group structure of GRBs. Somewhat later several authors [8, 9, 10, 11, 12] 
included more physical parameters into the analysis of the bursts (e.g. peak-fluxes, fluences, hardness ratios, etc.). A 
cluster analysis in this multidimensional parameter space suggests the existence of the third ("intermediate") group as 
well [7, 8, 13, 14, 12]. The physical existence of the third group is, however, still not convincingly proven. However, 
the celestial distribution of the third group is anisotropic [15, 16, 17]. All these results mean that the existence of the 
third intermediate group in the BATSE sample is acceptable, but its physical meaning, importance and origin is less 
clear than those of the other groups. Hence, it is worth to study new samples like the Swift data. 

In Sect. 2 we discuss the method used in the paper In Sect. 3 uni-, bi-, tri- and tetra-modal log-normal fits made 
by using the maximum likelihood method are discussed. In Sect. 4 one thousand Monte-Carlo simulations are shown 
concerning the significance of the fits. In Sect. 5 we discuss some further details. The conclusions are given in Sect. 6. 



THE METHOD 

There are several methods to test significance. For example the method which we used in our first paper [6] to 
analyze the Tgo distribution of the BATSE bursts is not useful here, because of the small population of short bursts in 
the Swift sample. 

In the Swift BAT Catalog [18] there are 237 GRBs, of which 222 have duration information. Fig. 1. shows the 
logTgo distribution. To use the method one has to bin the data. If the number of counts within some bins is small 
the method is not applicable. The Maximum Likelihood (ML) method is not sensitive to this problem, therefore for 
the (Swift) BAT bursts the maximum likelihood method is much more appropriate. 



The ML method assumes that the probability density function of an x observable variable is given in the form 
of g{x,pi, ...,pk) where pi,...,pk are parameters of unknown value. Having N observations of x one can define the 
likelihood function in the following form: 

N 

l = Y[si^i^Pl'--^Pk), (1) 

;=1 

or in logarithmic form (the logarithmic form is more convenient for calculations): 

N 

L = logl = Y,^og{g{xi,pu--,Pk)) (2) 

(=1 

The ML procedure maximizes L according to /? i ,...,/?<. . Since the logarithmic function is monotonic the logarithm 
reaches the maximum where I does as well. The confidence region of the estimated parameters is given by the following 
formula, where L^ax is the maximum value of the hkehhood function and Lq is the hkehhood function at the true value 
of the parameters [19]: 

2{L,nax-Lo)^Xk, (3) 



LOG-NORMAL FITS OF THE DURATION DISTRIBUTION 

Similarly to [20] we fit the logTgo distribution using ML with a superposition of k log-normal components, each of 
them having 2 unknown parameters to be fitted with A'^ = 222 measured points in our case. Our goal is to find the 
minimum value of k suitable to fit the observed distribution. Assuming a weighted superposition of k log-normal 
distributions one has to maximize the following likelihood function: 

Lk = ijlog w,/,(x,-,log7],(7;)j (4) 
where w/ is a weight, // a log-normal function with log 7} mean and CJ; standard deviation having the form of 

1 / (x-log7})2\ 
// = ^exp (5) 



and due to a normaUzation condition = A'^. We used a simple C++ code to find the maximum of Lj^. Assuming 
only one log-normal component the fit gives Li,^ = 951.666 but in the case of k=2 one gets Limwc = 983.317 with 
the solution displayed in Fig. 1. 

Based on Eq. (3) we can infer whether the addition of a further log-normal component is necessary to signifincantly 
improve the fit. We take the null hypothesis that we have already reached the the true value of k. Adding a new 
component, i.e. moving from to ^ -|- 1, the ML solution of L^max change to 7-(/S:+i)maA:' t>ut Lq remained the same. In 
the meantime we increased the number of parameters by 3 (w^+i, logTj^j^i and a(j^^i-^). Applying Eq. (3) to bothLj,^ 
and we get after subtraction 

2(-f'(fc+l)max ~ -^fcrnoA:) ~ Xl ■ (6) 

For A: = 1 L2max is greater thanLi,„a;i by more than 30, which gives for an extremely low probability of 5.88 x 10^'^. 
This means that the two log-normal fit is really a better approximation for the duration distribution of GRBs than one 
log-normal. 

Thirdly, a three-log-normal fit was made combining three fi, functions with eight parameters (three means, three 
standard deviations and two weights). The highest value of the logarithm of the likelihood (L^max) is 989.822. For 
two log-normal functions the maximum was L2max — 983.317. The maximum thus improved by 6.505. Twice this is 
13.01 which gives us the probability of 0.461% for the difference between L2max and L'^max is being only by chance. 
Therefore there is only a small chance the third log-normal is not needed. Thus, the three-log-normal fit (see Figure 
1.) is better and there is a 0.0046 probabihty that it was caused only by statistical fluctuation. 

One should also calculate the hkehhood for four log-normal functions. The best logarithm of the ML is 990.323. It 
is bigger by 0.501 than it was with three log-normal functions. This gives us a low significance (80.1%), therefore the 
fourth component is not needed. 
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FIGURE 1. Fit with two and three log-normal component for the duration distribution of BAT bursts. 



1000 MONTE-CARLO SIMULATIONS USING THE TWO-COMPONENT FIT 



We can check the 0.0046 probability, which we get for the maximum likeHhood calculation, using a Monte-Carlo (MC) 
simulation and adopting the following procedure. Take the two-log-normal distribution with the best fitted parameters 
of the observed data, and generate 222 numbers for Tgo whose distribution follows the two-log-normal distribution. 
Then find the best likelihood with five free parameters (two means, two dispersions and two weights; but the sum of 
the last two must be 222). Next we perform a fit with the three-log-normal distribution (eight free parameters, three 
means, three dispersions and two independent weights). Finally, we take the difference between the two logarithms of 
the maximum likelihoods that gave one number in our MC simulation. 

We have carried out this procedure for 1000 simulations each with 222 simulated log Tgos. There were 8 cases when 
the log-hkeUhood difference was more than the one obtained for the BAT data (6.505). Therefore the MC simulations 
confirm the result obtained by applying Eq. (6) and give a similar (0.8%) probability that a third group is merely a 
statistical fluctuation. 



DISCUSSION 

It is possible that the fit using three log-normal functions is accidental, and that there are only two types of GRBs. 
However, the probabiUty that the third component is only a statistical fluctuation is 0.5-0.8 %. 

One can compare the burst group weights with previous results. BAT sensitivity is different to BATSE sensitivity 
[21, 22]. BAT is more sensitive at low energies which means it can observe more X-ray flashes and soft bursts and 
probably fails to detect many hard bursts (typically short ones). Therefore one expects more long and intermediate 
bursts and fewer short GRBs. In the BAT data set there are only a few short bursts. Our analysis could only find 16 
short bursts (7%). The robustness of the ML method is demonstrated here because a group with only 7% weight is 
identified. Previously in the BATSE database intermediate bursts were identified by many research groups. However, 
in this class different frequencies were found representing 15-25 % of BATSE GRBs [7, 8, 13, 14, 23]. 



CONCLUSION 



1. Assuming that the Tgo distribution of the short and long GRBs is log-normal, the probabiUty that the third group 
is a chance occurance is about 0.5-0.8 %. 

2. Although the statistics indicate that a third component is present, the physical existence of the third group is still 
debatable. The sky distribution of the third component is anisotropic as proven by [15] and [16]. Alternatively 
[8] believe the third statistically proven subgroup is only a deviation caused by complicated instrumental effects, 
which can reduce the duration of some faint long bursts. This paper does not deal with this particular effect, 
however the previously studied BATSE sample shows a similar group structure. This agreement suggests that the 
third component is possibly real, not an instrumental effect (the BATSE detectors and the Swift BAT are different 
kinds of instruments). 

3. The observed frequencies in the three classes are different for BATSE and BAT. Both samples are dominated, 
however, by the long bursts. The short bursts are less populated in BAT than in BATSE but the intermediate group 
is more numerous. This is understandable, since BAT is less sensitive in high energy than BATSE was and more 
sensitive in low energy and short bursts are the hardest group and intermediate ones are the softest. Therefore 
BAT can observe more intermediate bursts and much fewer short ones than BATSE did. 

4. The existence and physical properties of the intermediate group need further discussion to elucidate the reality 
and properties of this class of GRBs. 
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